![]() Method for measuring the unbalance of elastic shaft rotors by means of displacement sensors
专利摘要:
In a method for determining unbalance compensation of resilient shaft rotors (1), the rotor (1) is rotatably housed in two bearing devices. A speed sensor (4) measures the speed of the rotor (1) and a radial movement of the rotor (1) is measured at different rotational speeds of the rotor by means of displacement sensors (3) at measuring points (6) during an unbalance measurement cycle or during several unbalance measurement cycles. The entered measurements are transmitted to an operating unit (5) which determines the eccentricity values corresponding to the measuring points (6) by means of an extension of the influence coefficient method, so that unbalances by plane and eccentricities per measuring point are determined for each cycle. Figure for the abstract: Fig. 1 公开号:FR3077638A1 申请号:FR1901097 申请日:2019-02-05 公开日:2019-08-09 发明作者:Kai (Dr.) Trukenmüller 申请人:Schenck RoTec GmbH; IPC主号:
专利说明:
Description Title of the invention: Method for measuring the unbalance of rotors with elastic shaft using displacement sensors The invention relates to a method for determining the unbalance of rotors with elastic shaft using displacement sensors at points of measurement whose eccentricity with respect to the bearing surface of the rotor is unknown. To balance rigid rotors, it is generally sufficient to apply correction masses in two axially separated planes. The correction masses, and therefore also the rotor unbalance, are independent of the operating speed. The balancing speed can be freely chosen. It is often chosen to be as low as possible and is often lower than when used later. On the other hand, in the case of rotors with elastic shaft, the concept of unbalance must be generalized. The number of planes required can be more than two and the axial position of the planes is of great importance. Whether a rotor can be considered rigid or an elastic shaft depends not only on the mechanical properties of the rotor, but also on its operating speed. The more the operating speed increases, the more important the elastic properties play. Document DE 41 33 787 A1 discloses a balancing process for elastic rotors which, without a test mass cycle, determines the correction masses necessary to balance the unbalance of the rigid body and the deflections of the rotor with elastic shaft. First, at least one unbalance measurement is determined during an unbalance measurement cycle at a speed at which the rotor exhibits behavior with a rigid body. At least one other unbalance measurement is determined per bearing plane and per mode suitable for correcting at a speed situated in the zone of the modal speed to be taken into account. The unbalance correction is calculated in an operating unit using the unbalance measurements obtained and characteristic data specific to the rotor or the bearing to compensate for the rigid body unbalance and the proper mode component to be taken into account, a force footprint, independent of speed and constant, being calculated for each bearing plane, which also describes the unbalance effect of the elastic rotor. In this process, it is necessary to enter the specific characteristic data of the rotor and of the bearing, such as the data of dimension, shape and material used of the rotor to be balanced. Document DE 40 19721 A1 describes a balancing method for an elastic rotor which operates at or near critical speeds, with a correction in three or more correction planes and the use of combinations of general unbalance distributions of the rotor and predetermined eigen modes, without actually balancing the rotor at critical speeds. During an unbalance measurement cycle at low rotational speed, corrections are usually determined which are made in two rotor correction planes. In addition, a third correction is made on the rotor which is proportional to the first and to the second correction and to the unbalance / form of mode combination. After another low-speed unbalance test cycle, which serves as the test mass cycle for the third correction, corrections are again determined for the first and second correction plan and used for the final balancing of the rotor. This balancing process does not make it possible to minimize the maximum admissible correction masses nor to minimize the admissible residual unbalance values. The unbalance measurement is often carried out indirectly by the vibration of the frame. This means that we measure the vibrations which are transmitted by the shaft to the fixed structure of the frame via bearings. The vibrations of the frame are generally measured by acceleration sensors (eg piezoelectric vibration sensors) or by vibration speed sensors. In particular for balancing installations at high speed or in the event of balancing in service of the rotors of power stations, such as turbines or generators, displacement sensors are also used which measure directly the radial movement of the shaft instead of the vibrations of the frame. The displacement sensors are for example capacitive or inductive displacement sensors (eddy current sensors) or laser triangulation sensors. In flexible bearings such as leaf or air bearings, the bearing gap is large relative to the eccentricity of the center of gravity of the rotor subjected to unbalance. The rotor rotates around its axis of inertia, which rotates eccentrically in the bearings. The forces transmitted to the surrounding structure (eg the frame) by the flexible bearings are relatively small. This has the consequence that it is certainly possible to measure with the aid of displacement sensors the axis of rotation extending eccentrically, but the vibrations transmitted to the structure are however too weak for an accurate determination of the unbalance by measurement. vibrations of the frame. The effect will be enhanced if the frame is relatively heavy compared to the rotor. This is particularly the case for units that have small rotors and a relatively large frame (for example with input and output lines or stator windings for an electric drive or generator). The influence coefficient method is a method of calculating the unbalance which is based on measurement cycles with known unbalances (eg reference work by W. Kellenberger; Elastisches Wuchten, Berlin 1987, pages 317-325). It is mainly used to measure the unbalance of rotors with elastic shaft. For this, it is assumed that there is a linear relationship between the applied unbalances and the oscillatory responses at the measurement points. The test unbalance per plane is noted as a complex number, where the amplitude and the angle of the complex number correspond respectively to the amplitude of the unbalance and to the angular position on the rotor. For each test cycle, the vibrations are measured if possible in several positions and at several speeds. The measurement can also be taken when the rotor is started, then assigned to fixed speeds or speed intervals. The type of measurement (e.g. acceleration at bearing positions or measurement of the radial position of the shaft using displacement sensors) is not important for the formulation of the influence coefficient method. From the time series of sensor values, the first order is evaluated in relation to the current rotation frequency. This requires a speed or angle sensor that measures a mark on the rotor. A complex number is also formed from the signal amplitude and the phase of the vibration sensors relative to the speed sensor. The relationship between the unbalance and the vibration signal can then be described using a system of linear equations with complex numbers: [0007] [Math.l] S - [0008] Here, 5 is the vector of the complex amplitudes of the signal. It contains the values of all the sensors and speeds used. It therefore has the length N = (number of speed support points) x (number of sensors). The unbalances are gathered by plane of unbalance in the vector of unbalances u. For P unbalance planes, u has the length P. The matrix of the influence coefficients K has the dimension N x P. Both 5 and K and u are complex quantities. The influence coefficients method can be solved, for example, by determining the matrix K. This is generally done by P + l test cycles, where a known single unbalance is applied to each cycle in a plane of different imbalance. In addition, it is possible to carry out a cycle without test mass, a cycle which is called the initial cycle (Urlauf). Then the matrix K of the influence coefficients is determined by subtracting the initial cycle from all the other test cycles. However, this is only one of the possible methods for determining K. When K is known, it is possible to calculate the unbalance in all planes for any cycle of the rotor. This is done by solving the system of linear equations according to the unknown unbalance vector u. In general, there are more measurement points N than unknown unbalances P, so that the system of equations is overdetermined. In this case, it can be solved using the linear correction calculation (least squares method). However, there are other methods which solve the system using optimization methods and, where appropriate, additional boundary conditions. Once the system has been resolved, the imbalances u are generally compensated in each plane, for example by the application of counterbalances or by removal of material. The vibration can then be measured again and, if necessary, the rotor can be balanced in additional steps. If displacement sensors are used in the influence coefficient method, the objective is then to reduce the bending of the shaft caused by the unbalance. This becomes problematic if the measuring point on the rotor, where the bending is measured, already has eccentricity without rotation, for example due to a bent shaft or an impact on the measuring surface. Significant errors can be generated if eccentricity errors are ignored and the influence coefficient method is applied according to the known method. In this case, the constant measurement error is interpreted as the vibratory displacement of the shaft which must be compensated for by balancing. The objective of the invention is to provide a method of measuring unbalance in which an eccentricity in at least one measurement point is taken into account. This object is achieved by a method for determining an unbalance of rotors with elastic shaft by means of displacement sensors at measurement points which have an unknown eccentricity with respect to the bearing surfaces, method in which the rotor is housed. rotatable in bearings, a speed sensor measures the rotor speed, a radial movement of the rotor is measured by means of displacement sensors at measuring points during an unbalance measurement cycle or during several unbalance measurement cycles at different rotational speeds of the rotor, the measurements entered are transmitted to an operating unit, the operating unit using the measurements entered to solve a system of equations which, in addition to the unbalance, includes the unknown eccentricity, so that the imbalance per plane and the eccentricity per measurement point are calculated. Using the method according to the invention, it is possible to determine the eccentricity of the measurement point during the unbalance measurement cycle proper and to take it into account when calculating the unbalance. This means that no separate measurement cycle is necessary to determine a possible eccentricity at low speed. The unbalance values corresponding to the measurement points can, for example, be determined by extending the method of the influence coefficients. However, other optimization methods can also be used. If the measurement points have unknown eccentricities, these are contained in the complex amplitudes of the signal in the form of errors. The system of equations of the method of influence coefficients can be advantageously extended by [Math.2] s = K ii 4- E e [0017] where e is the vector of the unknown eccentricities. The length of the vector is equal to the number of measurement points. The matrix E consists only of zeros and ones. There is only one in each column. It is located in column number k, which corresponds to the measurement point of this line: with 1 if signal s t from the sensor at the position, otherwise 0. lo, [0018] The extension of the influence coefficient method can be formulated as follows: [Math.3] s - 4 λ [0020] For this purpose, the two KetE matrices can be written next to each other as follows: [Math.4] A = [ALE] at A. The vector x can include unknown imbalances and eccentricities. [Math.5] The determination of the matrix of the influence coefficients K can be carried out according to the known method of the influence coefficients. As differences are always formed, the basic eccentricity has no influence. In one embodiment, the extended influence coefficient method can be used to calculate the unbalance for any measurement cycle. It is preferable that the unbalance is measured in a speed range in which the rotor has elastic deviations. When balancing at low speed, it is not always possible to rotate the rotor around the axis that the component will occupy in the subsequent mounting of the rotor. This can be disadvantageous during subsequent operation of the rotor, since even small eccentricity errors are sufficient to cause significant unbalance values on the mounted rotor. In a preferred configuration, balancing takes place in a speed range in which the rotor already has an elastic shaft behavior. This simulates the normal operation of a rotor system rotating at high speed and balances the rotor at speeds up to service speed. To measure the unbalance, it may be preferable to accelerate the rotor to predefined measurement speeds. This makes it possible to accurately determine the effects of a rotor unbalance and, for example, a test unbalance at a desired speed. It may be advantageous if the measurements are made during acceleration or deceleration of the rotor, so that it is not necessary to carry out other test cycles and that the balancing can be done effectively. It may also be preferable to determine the unbalance and the eccentricity in the service speed range of the rotor. In this case, the rotor is accelerated to a speed which advantageously lies at its service speed or at least close to its service speed. With elastic shaft rotors precisely, balancing at low speed is not satisfactory. On the other hand, with the preferred embodiment of the method, the rotor is accelerated at speeds within its operating speed range, so that the unbalance is measured in a speed range which reflects the actual operation of the rotor and which allows a reliable balancing. In a preferred embodiment, the rotor is housed in flexible bearings. Rigid bearings are often used to transmit unbalance vibrations to the frame and to measure them. For certain rotor systems, it is however advantageous, even necessary, to use flexible bearings. Thanks to the preferred embodiment of the method, in which the rotor is housed in flexible bearings and the measurement is carried out using displacement sensors, the previous measurement methods can be considerably improved and the eccentricity can be determined by more unbalance. It may be advantageous to measure vibrations at other measurement points by means of a sensor element. In the case of vibration transmission, vibration sensors can be used in addition to displacement sensors, the measurements of which can be used to solve the system of equations. The invention is explained in more detail below using an embodiment of the invention shown in the figure. Figure [fig.l] shows a rotor in an unbalance measuring device. Figure 1 schematically shows an unbalance measuring device in which a driven rotor (1) has bearing surfaces (2) and is housed, for example, in two bearings. The rotor (1) can also be housed in more than two bearings. The bearings can be flexible bearings, such as leaf bearings or air bearings. Two displacement sensors (3) measure radial movements of the rotor (1), that is to say its bending. A speed and angle sensor (4) measures a mark on the rotor (1). The measurements captured by the sensors (3, 4) are transmitted to an operating unit (5). Displacement sensors (3) are used which determine an eccentricity of at least one measurement point (6), or of a measurement surface, which has an unknown eccentricity with respect to the bearing surfaces (2 ). Due to the soft bearings, no significant vibration occurs which could be used to calculate the unbalance. If necessary, additional sensors can be used to measure for example a vibration of the frame or a vibration of the foundations, which in turn can be used to calculate the unbalance. The present method can be used advantageously when the measurement of the eccentricity cannot be carried out at the positions of the bearings because the bearings are closed or inaccessible. It is also possible to take the measurements directly next to the bearings, at a point which has only a slight eccentricity with respect to the bearing. However, this is only possible if the bearing positions are accessible to the displacement sensors (3). However, the rotor (1) is generally designed to be as compact as possible. Therefore, it may not be possible to place sensors at bearing positions, particularly in the case of small units, due to spatial constraints. However, it is possible to easily measure the imbalance of compact rotor systems thanks to the preferred method and the placement of the displacement sensors (3) on measurement points whose radial eccentricity is unknown relative to the bearing surfaces (2) . A drive device for the unbalance measurement device accelerates the rotor (1) at appropriate speeds, the drive device can also be controlled so that it is possible to preset fixed balancing speeds . The measurements to be entered can also be entered in normal operation, so that, for example, when accelerating the rotor (1) at balancing speeds which can be preset, the corresponding measurements are entered and transmitted to the unit. operating (5). The rotor (1) can also be accelerated at speeds which lie within its operating speed range. It is possible with the method of balancing rotors (1) having any number of balancing planes. Often, balancing plans are imposed by design, as in turbochargers or turbine rotors, and allow only a maximum allowable balancing mass for reasons of size or strength. The rotor (1) is accelerated to a preset rotational speed n at which it exhibits the behavior of a rotor with an elastic shaft. When the preset speed n is reached, the sensors (3, 4) take measurements and transmit them to the operating unit (5). The measurement can also be carried out at several rotational speeds n. In the case where the eccentricity error of the measurement points (6) relative to the bearing surfaces (2) is known, it can be taken into account mathematically in the measurement signals and the known method of influence coefficients can, for example, be used to calculate the unbalance. However, it is not always possible to measure eccentricity. The bearings are generally no longer accessible in the assembled state. When the measurement is made before mounting, it is not possible to guarantee that eccentricity errors will be retained during mounting. Often, the eccentricity of the measurement surface (6) is determined at low speed and the measured values are subtracted from the measurements at high speed (so-called runout compensation / runout compensation). However, this only works with rigid bearings, such as roller bearings. In the case of hydrodynamic plain bearings, different axes of rotation (orbits) often occur at different speeds of rotation, so that an eccentricity error cannot be unequivocally defined. Other bearings, for example leaf bearings, need a minimum speed to be able to operate stably. In such cases, it is often not possible to slowly rotate the rotor. It is accelerated very quickly to the minimum speed of rotation which is beyond the natural frequencies of the rigid body of the rotor-bearing system. Using the method according to the invention, the unbalance measurement can be easily carried out taking into account the eccentricity of such rotor systems. An unbalance of elastic shaft rotors (1) can be determined by means of displacement sensors (3) on measuring points (6) having an unknown eccentricity with respect to the bearing surfaces (2), the rotor (1) being rotatably housed in bearings, the speed sensor (4) measuring the speed of the rotor (1), a radial movement of the rotor (1) being measured by means of displacement sensors (3) at measuring points (6) during a measurement cycle or several unbalance measurement cycles at different rotor speeds, and the captured measurements being transmitted to an operating unit (5), the operating unit (5) using the entered measurements to solve a extension of the influence coefficients method according to the system of equations, with vector 5 which includes measures and velocities, with =. [K E. with the matrix of influence coefficients K of dimension, v .... "and the matrix r Ί and with 1 if J lo, signal from the sensor at the position, otherwise 0, with the vector x including the unbalances and eccentricities unknown with , so that for each measurement cycle it is possible to determine unbalances per plane and eccentricities per measurement point. A central idea of the process is that the eccentricity, in addition to the unbalance to be determined, is contained as an unknown in the system of equations, the system of equations being able to be solved by optimization methods. The use of an extended coefficients of influence method is only a preferred embodiment. The method has proved particularly advantageous for rotor systems rotating at high speed in small units such as electric compressors, vacuum cleaner motors, pumps or micro gas turbines. In these systems, the rotor is often housed in flexible bearings, such as air, gas, foil or magnetic bearings. With the method according to the invention, it is possible to measure the unbalance and to determine the eccentricity, the displacement sensors advantageously measuring the bending of the rotor shaft outside the bearing positions.
权利要求:
Claims (1) [1" id="c-fr-0001] claims Method for determining an unbalance of elastic shaft rotors (1) by means of displacement sensors (3) at measurement points (6) whose eccentricity with respect to the bearing surface (2) is unknown, in which the rotor (1) is rotatably housed in bearings, a speed sensor (4) measures the speed of the rotor (1), a radial movement of the rotor (1) is measured by means of displacement sensors (3) measurement (6) during an unbalance measurement cycle or during several unbalance measurement cycles at different rotational speeds of the rotor, the captured measurements are transmitted to an operating unit (5), the operating unit using the measurements entered to solve a system of equations which, in addition to the unbalance, includes the unknown eccentricity as unknown, so that the unbalance by plane and the eccentricity by measurement point are calculated. Method according to claim 1, characterized in that the measurement is carried out in a speed range in which the rotor (1) has elastic deviations. Method according to claim 1, characterized in that the rotor (1) is accelerated at predefined measurement speeds. Method according to claim 1, characterized in that the measurements are taken when the rotor (1) is accelerated or decelerated. Method according to claim 1, characterized in that the unbalance and the eccentricity are determined in the operating speed range of the rotor (DProcess according to one of the preceding claims, characterized in that the rotor (1) is housed in flexible bearings. Method according to one of the preceding claims, characterized in that vibrations are measured at other measuring points by means of a sensor element.
类似技术:
公开号 | 公开日 | 专利标题 FR3077638A1|2019-08-09|Method for measuring the unbalance of elastic shaft rotors by means of displacement sensors EP1382858B1|2008-01-09|Distribution method for rotor blades of a turbomachine EP1439010B1|2018-11-14|Unbalanced dynamic vibration generator EP0505976B1|1994-09-21|Device to compensate for the vibratory force or couple felt by a body FR2548369A1|1985-01-04|COMPENSATED RHEOMETER EP2936056B1|2017-08-09|Gyroscope with simplified calibration and simplified calibration method for a gyroscope EP2992298B1|2020-03-18|Mems balanced inertial angular sensor and method for balancing such a sensor FR2580401A1|1986-10-17|APPARATUS AND METHOD FOR MEASURING THE VISCO-ELASTIC PROPERTIES OF MATERIALS FR2902879A1|2007-12-28|ORTHOGONAL RHEOMETER EP0337040A1|1989-10-18|Device for compensating a vibrational force or a vibrational torque created by a body WO2010142866A1|2010-12-16|Dynamic load bench FR2844877A1|2004-03-26|Turbocompressor rotor balancing method in which both primary and secondary resonance effects are balanced simultaneously EP1355139B1|2013-12-04|method and device for correcting an unbalance EP0116253B1|1987-09-02|Apparatus for measuring the amplitude and position of an unbalance of a turning system EP1371960A3|2004-05-12|Dynamic balance testing machine FR2677121A1|1992-12-04|Method and device for balancing a flexible rotor EP1326065A1|2003-07-09|A method of measuring the unbalance of rotors EP2558819B1|2014-07-16|Gyroscopic measurement in a navigation system NL8300126A|1983-08-01|BALANCING A ROTATION BODY. EP0702420A1|1996-03-20|Surface acoustic wave motor EP3287237B1|2019-08-14|Screwing device with optimised output torque measurement, and method for determining the corresponding output torque EP0336795B1|1991-09-11|Mechanical type composite inertial sensor EP0810418A1|1997-12-03|Apparatus for measuring rotation FR2825485A1|2002-12-06|METHOD AND DEVICE FOR CONTROLLING THE ANGULAR SPEED OF A LITTLE CUSHIONED ELECTROMECHANICAL CHAIN EP2406878B1|2017-05-03|Device for quickly generating a torque for an extended dynamic range with low inertia
同族专利:
公开号 | 公开日 CN110118632A|2019-08-13| DE102018102751B3|2019-02-21| US20190242774A1|2019-08-08| FR3077638B1|2021-07-23| US10823632B2|2020-11-03|
引用文献:
公开号 | 申请日 | 公开日 | 申请人 | 专利标题 US3096639A|1960-02-24|1963-07-09|United Aircraft Corp|Method and apparatus for compensating for eccentricity| AT245825B|1960-09-13|1966-03-25|Schenck Gmbh Carl|Machine for determining the unbalance in balancing bodies| US4135244A|1977-07-20|1979-01-16|Spectral Dynamics Corporation|Method and apparatus for programming and reducing runout in rotating shaft vibration signal detection| US4170896A|1978-03-20|1979-10-16|Summa Corporation|Balancing of high-speed shafts| FR2506455B1|1981-05-21|1984-06-01|Elf Aquitaine| FR2538903B1|1983-01-03|1985-08-02|Snecma|APPARATUS FOR MEASURING THE AMPLITUDE AND THE ANGULAR POSITION OF A LOOP OF A ROTATING SYSTEM| US5214585A|1989-06-30|1993-05-25|General Electric Company|Balancing method and product| US5140534A|1990-07-20|1992-08-18|Westinghouse Electric Corp.|Centerless runout and profile inspection system and method| DE4133787C2|1991-10-11|2002-06-20|Schenck Rotec Gmbh|Balancing method for the determination of the balancing masses for an elastic rotor on a force-measuring balancing machine and device for carrying out the method| US6757636B2|2002-07-19|2004-06-29|Alstom Technology Ltd.|Computerized electronic runout| GB2398882B|2003-02-27|2007-09-05|Rolls Royce Plc|Rotor balancing| CN1632495A|2004-12-08|2005-06-29|天津大学|Flexible rotor dynamic balance determination apparatus based on transfer function| US7500432B2|2005-10-28|2009-03-10|Van Denend Mark E|Apparatus and method for balancing a printing roller having an image producing area on its outer surface| GB2438642B|2006-06-01|2008-04-09|Rolls Royce Plc|Rotor unbalance correction| CN101625277B|2008-07-07|2011-07-27|西门子公司|Method and device for quantitatively detecting nonequilibrium state and method for detecting clamping state of workpiece| CZ303892B6|2011-09-21|2013-06-12|Doosan Skoda Power S.R.O.|Method of determining current eccentricity of rotating rotor and rotating rotor eccentricity diagnostics| CN103076163B|2011-12-06|2016-02-24|西安交通大学|A kind of on-line testing method of Rotor-Bearing System characterisitic parameter| CN103776587B|2014-01-28|2016-03-16|郭卫建|Determine the method for the amount of unbalance of rotor| CN104501714B|2014-12-29|2017-03-15|贵州电力试验研究院|Turbine rotor bias azimuth on-Line Monitor Device and its monitoring method| CN105021349A|2015-05-19|2015-11-04|郭卫建|Method for obtaining unbalance of rotor| US10753817B2|2017-03-06|2020-08-25|Yaskawa America, Inc.|Testing apparatus, computer readable medium, and method for minimizing runout| CN107133387B|2017-04-10|2019-10-18|浙江大学|The imbalance compensation method of rotor unbalance coefficient variable step polygon iterated search|DE102018115363B3|2018-06-26|2019-11-07|Schenck Rotec Gmbh|Method for determining an imbalance of a wave-elastic rotor based on the deflection| CN111413030B|2019-01-07|2021-10-29|哈尔滨工业大学|Large-scale high-speed rotation equipment measurement and neural network learning regulation and control method and device based on rigidity vector space projection maximization| CN111413031B|2019-01-07|2021-11-09|哈尔滨工业大学|Deep learning regulation and assembly method and device for large-scale high-speed rotation equipment based on dynamic vibration response characteristics|
法律状态:
2020-02-20| PLFP| Fee payment|Year of fee payment: 2 | 2020-11-27| PLSC| Publication of the preliminary search report|Effective date: 20201127 | 2021-02-17| PLFP| Fee payment|Year of fee payment: 3 | 2022-02-21| PLFP| Fee payment|Year of fee payment: 4 |
优先权:
[返回顶部]
申请号 | 申请日 | 专利标题 DE102018102751.6A|DE102018102751B3|2018-02-07|2018-02-07|Method for measuring the imbalance of wave-elastic rotors by means of distance measuring sensors| DE102018102751.6|2018-02-07| 相关专利
Sulfonates, polymers, resist compositions and patterning process
Washing machine
Washing machine
Device for fixture finishing and tension adjusting of membrane
Structure for Equipping Band in a Plane Cathode Ray Tube
Process for preparation of 7 alpha-carboxyl 9, 11-epoxy steroids and intermediates useful therein an
国家/地区
|